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The Rashba spin-orbit coupling arising from structure inversion asymmetry couples spin and 
momentum degrees of freedom providing a suitable (and very intensively investigated) environment 
for spintronic effects and devices. Here we show that in the presence of strong disorder, non¬ 
homogeneity in the spin-orbit coupling gives rise to a finite spin Hall conductivity in contrast with 
the corresponding case of a homogeneous linear spin-orbit coupling. In particular, we examine the 
inhomogeneity arising from a striped structure for a two-dimensional electron gas, affecting both 
density and Rashba spin-orbit coupling. We suggest that this situation can be realized at oxide 
interfaces with periodic top gating. 

PACS numbers: 72.25.-b, 75.76.-l-j, 72.25.Rb, 72.15.Gd 


The spin Hall effect (SHE) [T] is the generation of a 
transverse spin current by an applied electric field with 
the current spin polarization being perpendicular to both 
the field and the current flow. Since the SHE allows the 
control of the spin degrees of freedom even without exter¬ 
nal magnetic fields (see e.g. Refs. [HE]), it has become 
a central topic in present spintronics research. [H [5] The 
microscopic origin of the SHE lies in the spin-orbit cou¬ 
pling (SOG), which in solid-state systems may be due to 
the potential of the ionic cores of the host lattice, the po¬ 
tential of the impurities and the confinement potential of 
the device structure. In a two-dimensional electron gas 
(2DEG), Bychkov and Rashba [6] have proposed that the 
lack of inversion symmetry along the direction perpen¬ 
dicular to the gas plane leads to a momentum-dependent 
spin splitting usually described by the so-called Rashba 
Hamiltonian 

H=- -harxz-p, (1) 

2m 

where p is the momentum operator for the motion along 
the plane hosting the 2DEG, say the xy plane, z is a 
unit vector perpendicular to it, t = is the 

vector of the Pauli matrices, and a is a coupling con¬ 
stant whose strength depends on both the SOG of the 
material and the field responsible for the parity break¬ 
ing. The Hamiltonian Q, which has been extensively 
used in the study of the 2DEG in semiconducting sys¬ 
tems, has been recently applied also to interface states 
between different metals [7] and between two insulating 
oxides [9HI2]. In the latter systems, higher mobilities, 
carrier concentration and SOC strengths have led to the 
expectation of observing stronger SOG-induced effects. 
The Hamiltonian 0 is deceptively simple, as one realizes 
when considering transport phenomena. In particular, 
the intrinsic universal SHE proposed in Ref. m turned 
out to be a non stationary effect, while under station¬ 
ary conditions cancellations occur, leading to a vanish¬ 
ing spin Hall conductivity (SHG) i.e., the coefficient 


relating the z-spin current in the y direction to the ap¬ 
plied electric field, Jy = Here, we show that 

this is only true for a spatially homogeneous a, while 
considering a space-dependent a{x, y) opens the way to 
a substantial SHE under stationary conditions, even in 
the presence of strong disorder. We shall first discuss 
from a general perspective how this comes about, and 
shall afterwards demonstrate numerically the effect in the 
presence of a spatially modulated SOG as it could be re¬ 
alized in the 2DEG at the interface of a LaAlOa/SrTiOa 
(LAO/STO) heterostructure, schematically depicted in 
Eig.[T] It is experimentally established that the Rashba 
SOG increases when the electron density in the 2DEG 
of these heterostructures is increased [HHiniiii]- Since 
the local electric field determining a is tightly related to 
the electron density m, one can naturally infer that the 
structure in Eig.[^ produces a modulation of the Rashba 
SOG. 
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FIG. 1. Schematic view of a possible device in which the 
SHE is enforced in the 2DEG at the interface of a LAO/STO 
heterostructure. The yellow stripes represent top-gating elec¬ 
trodes of width w and interspacing d. 

— General arguments — The interplay of the intrinsic 
Rashba SOG with the scattering from impurities makes 
the dynamics of charge and spin degrees of freedom in¬ 
trinsically coupled in a subtle way. This is especially evi- 
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dent in the vanishing of the SHC with homogeneity in the 
SOC. Notice that the spin current is a tensor quantity de¬ 
pending on the flow direction (lower index) and spin po¬ 
larization axis (upper index) and hence the spin current 
and the electric field are related by a tensor of third rank 
afj. For fixed polarization a = z, Onsager’s relations 
require the antisymmetry property cr^^, = — 

The vanishing of the SHC manifests via an exact 
compensation of the contribution originally proposed by 
Sinova et al. m by a further contribution, which arises 
by the coupling between the spin current and the spin 
polarization, which is induced in the plane perpendic¬ 
ularly to the applied electric field. This latter effect 
was almost simultaneously proposed by Edelstein [TS] 
and by Aronov and Lyanda-Geller [TB]. It consists of a 
non-equilibrium spin polarization due to the electric field 
Sq = —eaNorEx for a field along the x axis. Here e 
is the unit charge (—e for electrons), Nq = is the 

density of states per spin of the 2DEG described by the 
first term of Eq. Q and r is the elastic relaxation time 
introduced by impurity scattering. In such a 2DEG the 
standard Drude formula can be written via the Einstein 
relation as cr = 2e^NQD, with the diffusion coefficient 
D = VpTf2, vp being the Fermi velocity related to the 
Fermi energy Ep = mvp/2. To understand the origin of 
the compensation mentioned above, it is useful to start 
from a property of the Hamiltonian Q first pointed out 
by Dimitrova m, which relates the time derivative of 
the spin polarization to the spin current 

dtSy = -2mah-^J^. ( 2 ) 

Notice that such a relation is not changed by disorder 
scattering as long as the latter is spin independent. Dis¬ 
order is necessary to guarantee a stationary state which 
implies the left-hand side of Eq. (i) to be zero. Obvi¬ 
ously, the corresponding vanishing of the right-hand side 
entails a vanishing SHG when a is a constant. The spe¬ 
cific way in which the vanishing of the SHG occurs in 
a disordered 2DEG via the so-called vertex corrections 
[IHHSI] can be heuristically understood by describing the 
coupling between spin current and spin polarization as a 
generalized diffusion in spin space. By dimensional ar¬ 
guments spin current and spin density must be related 
by the factor Lso/tdp, where Lso = h/{2ma) is the 
spin-orbit precession length originating from the differ¬ 
ence of the Fermi momenta of the two branches of the 
spectrum of Eq. 0. while Tp)p is the Dyakonov-Perel 
spin relaxation time due to the interplay of SOG and dis¬ 
order scattering. In a disordered 2DEG, Tpp is related 
to Lso by the diffusion coefficient L’gQ = Dtdp, so that 
one obtains for the spin current 

= 2mah-^DSy + af E^. (3) 

In the above equation, which can be rigorously derived 
[53] , the quantity is the intrinsic contribution of 


Ref. [T3| in the diffusive regime appr/h <C 1, whereas 
the term proportional to corresponds to the vertex- 
correction contribution mentioned above. Given the ex¬ 
pression for (Tg^ = (e/87r)(2T/r£)p) derived in Ref.[T31it 
is now apparent that, if we replace with the Edelstein 
result, the spin current in Eq. ([^ vanishes, consistently 
with the stationarity requirement derived from Eq. (|^. 
The key observation is that this compensation does not 
necessarily occur for an inhomogeneous SOG, where it is 
no longer possible to express the time derivative of the 
spin polarization in terms of the spin current. In such a 
situation a spin current becomes possible under station¬ 
ary conditions. 

In order to illustrate the physical mechanism by which 
an inhomogeneous Rashba SOG leads to a finite SHE it 
is useful to consider a single-interface problem, which is 
described by Eq. Q with the replacement a —>■ a(x) = 
0(x)ap+0(—x)a- (with a+ > q;_). Glearly as x —)■ ± 00 , 
one recovers the uniform case with complete cancellation 
of the spin current for the Rashba model with couplings 
a±. On both sides of the interface, the y-spin polariza¬ 
tion obeys a diffusion-like equation with L±, the corre¬ 
sponding spin-orbit lengths in the two regions. One can 
then seek a solution of the form 

Sy{x)=9ix) (5o.+ +<5s+e-"/^+) 

+ 0i-x) (5o._+(5s_e"/^+) , 

where Sq ± are the asymptotic values of the y-spin po¬ 
larization at ±(X). The constants must be deter¬ 
mined by imposing the appropriate boundary conditions 
at X = 0. As a result, the spin current is exponentially 
localized near the interface, where the spin polarization 
S'^(x) must interpolate between the two asymptotic val¬ 
ues and there is no longer complete cancellation between 
the two terms of Eq. 

One can imagine to generalize this analysis to a series 
of interfaces apt to describe a periodic modulation of the 
SOG. Expectedly, if the spin-orbit length is larger than 
the distance between two successive interfaces, the spin 
Hall current should become practically uniform. 

— The model and its numerical solution — The previ¬ 
ous arguments within the diffusive limit are now substan¬ 
tiated by numerical results for a microscopic 2D lattice 
model (size x Ny) with inhomogeneous Rashba SOG 
described by the Hamiltonian 

H = J2 *L4aCja + - p)cl^Cia + ( 4 ) 

ij(T i(7 

where c|^(cicr) creates (annihilates) an electron with spin 
projection a on the site identified by the lattice vector i?i, 
the first term describes the kinetic energy of electrons on 
a square lattice (lattice constant a, only nearest-neighbor 
hopping: = —t for \Ri — Rj\ = a) and in the second 
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term Vi is a local disorder potential with a flat distribu¬ 
tion —Vo <Vi<Vo, and ^ is the chemical potential. The 
last term is the lattice Rashba SOC, 


ttRSO _ _• \ ^ 


iaa' 


^i+x,cr'’'a' 




i(7(7' 


‘ a a' 


J 


and the coupling constants ai^i+^/y > 0 are now defined 
on the bonds. Note that for constant a = Ui^i+x/y Eq- 
@ takes the usual form [2nd term in Eq. Q] in momen¬ 
tum space [23]. 

One can show m that a “continuity equation” for the 
local spin density Sf, (a dot stands for time derivative) 

~f [div J^ji -|- -|- o'i—y^iJi_y i = 0, (5) 


holds. For a homogeneous Rashba SOC, where 
[divJ^'ji = 0, Eq. ([^ corresponds to Eq. ([^ and im¬ 
plies that the total x-spin current has to vanish under 
stationary conditions. On the contrary, when a varies in 
space, a cancellation occurs [30] between div and the 
last two terms of Eq. ([^, so that the stationarity condi¬ 
tion S = Q can be fulfilled without the vanishing of . 
This is also clear for a system with periodic boundary 
conditions, where the total divergence of any current has 
to vanish, i.e.. 


^ ~ ^ ' {^^i,i+vJi,i+y + ^i-y,iJi-y,i^ ■ ( 6 ) 

i i 

Clearly, the left-hand side can vanish without implying 
= 0, because, if a is inhomogeneous, Eq. (H can 
be fulfilled with alternating signs of in regions with 
different a [see the top panel of Fig. 4 in Ref. EOj. 

We exemplify the situation for an inhomogeneous 
Rashba SOC which varies along the x direction form¬ 
ing a superlattice with d = 20a and w = 10a [see Fig.[^. 
In the regions of width d, a = ao is smaller than in the 
regions of width w, where a = ai > ao- The inhomogene¬ 
ity in ai^i+x,y leads to a concomitant charge modulation 
which is shown in Fig. ia). We have diagonalized the 
Hamiltonian Q and calculated the SHC from the Kubo 
formula (see, e.g., Ref. [T3|l tr®^ = cr®j^ with 

E-m 

Here we have taken the limit of zero temperature and 
r; —>■ 0 is a small regularization term which can be inter¬ 
preted as an inverse electric-field turn-on time [27] . Note 
that for the striped system one has already spin currents 
Jt^+y flowing in the ground state [2S] (see, e.g., Fig. 3 in 
Ref. liiOjl whereas the electrically induced spin current is 
Ji'i+y = Ej] f^tfEx- Thus Eq. ^ can be split into a 



FIG. 2. Main panel: Spin Hall coefficient and “station¬ 
arity” parameter 7 [both in units of l/(87r)] as a function of 
fj, for a homogeneous system with a = 0.51 (black solid fine) 
and a striped system {a‘^■. red solid and 7 : blue dashed) with 
modulated Rashba SOC as shown in the inset. Here the (red) 
dashed line displays the variation of a{x) = Ui^i+x = cei,i+y 
along the x direction for stripes along the y direction and 
width u) = 10 a separated by a distance d — 20 a. The Rashba 
SOC on the stripes is a(x) = ai — 0.81 while between the 
stripes a(x) = ao = 0.21. The black solid line (dots) in the in¬ 
set reports the charge profile at chemical potential fj, = —4.3 1. 
System size: 3060 x 3060 sites. 
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FIG. 3. a): (thick lines) and 7 (thin lines) as a function 

of density n (average number of electrons per lattice site) for 
a striped system and parameters are indicated in the panel, 
b): (thick fines) and 7 (thin lines) as a function of density 

but now with an additional modulation of the local chemical 
potential which is set to —0.5t on the stripes (ao = 0.21, 
ai = 0.8f). 


ground-state contribution (for which of course = 0 ) 
and a linear-response part. For the latter we define the 
quantity 7 = 2 which therefore also de¬ 

scribes the linear response of Ei Ef = — fEx to the ap¬ 
plied electric field and which in the following will be used 
to quantify the “stationarity” of the solution. In fact. 
Fig. demonstrates that for a constant a = 0.51 (black 
solid line) cr®'^ coincides with 7 and therefore the finite 
asH is a non-stationary result. On the other hand, the 
same panel also reports the results for the case Gq = 0 . 2 1 
and ai = 0.8 1. In this case one can see that for a non- 
negligible range of chemical potential near the bottom of 
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the band a substantial (red solid curve) is present 
while 7 = 0 (blue dashed curve), marking the occurrence 
of a SHE in stationary conditions. This indicates the rel¬ 
evant role of those states that are still extended along the 
y direction, while they are nearly localized inside the po¬ 
tential wells arising from the modulation of a. As shown 
in Fig.j^a), this situation occurs for increasingly large 
density ranges by increasing the inhomogeneity of a. We 
have also checked that the stationarity is not only global 
but that in the low density regime 5^ « 0 is fulfilled at 
each lattice site. 

We also investigated the effect of an inhomogeneous 
chemical potential as it is induced by the striped gating 
of Fig. [2 In particular, we shift the chemical potential 
downwards by = 0.5 1 on the sites below the gate (the 
regions of width w with a = oi = 0.81. From Fig.j^b) 
one sees that, although is reduced, it still remains 
substantial and the density range with a stationary SHE 
is even extended (black dashed curves). On the contrary 
one can see m that, in the absence of an inhomogeneous 
Rashba SOC, a simple charge modulation does not pro¬ 
duce any stationary SHE. 



FIG. 4. Main panel (a): Spin Hall conductivity as a function 
of chemical potential p and various values of disorder. Panel 
(b): as a function of p including error bars for Volt — 0.5 

and two values of y/t = 0 (circles) and y/t = 0.02 (triangles). 
Panel (c) reports the behavior of the “stationarity” parameter 
7 . Results in panels (a,c) are obtained for y — 0. 

Finally, we address the quite important issue of the 
robustness of SHE in the presence of disorder. Previ¬ 
ous analyses [261 El] showed that the SHC for a linear 
Rashba SOC is rapidly destroyed by disorder. This is 
easily understood because in the homogeneous Rashba 
SOC the SHE is a non-stationary effect which cannot 
survive the relaxation effect of disorder scattering. Here, 
instead, when 7 = 0 the SHE is present in a stationary 
state and disorder is much less effective in spoiling it. 
In the presence of a random potential of finite width Vq 
the calculations can only be carried out on smaller lat¬ 
tices (40 X 40 sites) where we consider stripes of width 


w = 4a and distance d = 4a (see Fig. [^. We follow the 
procedure described in Ref.|26]and diagonalize the Hamil¬ 
tonian for different disorder configurations and different 
twisted boundary conditions. For each concentration we 
consider 250 random boundary phases and 50 disorder 
configurations. 

As a striking result (main panel of Fig. we find 
that the average SHC at low densities (p < —4t) is not 
affected by disorder and only gets suppressed when the 
chemical potential is within the range of band states ex¬ 
tended both along x and y directions. In contrast, and as 
mentioned above, cr®^ vanishes for a homogeneous, lin¬ 
ear Rashba SOC [551 HZ] in the presence of disorder and 
for ?7 —>• 0 . 

It has been pointed out in Ref. m that the evalua¬ 
tion of CT®^ on finite lattices and taking the limit y ^ 0 
is complicated by strong fluctuations. These strong vari¬ 
ances in the SHC are exemplified in panel (b) of Fig. 
for Vo/t = 0.5 where we also show the corresponding re¬ 
sult for y/t = 0.02. The SHC in the low density regime 
is not dependent on the small y value but one observes a 
large reduction in the variance which becomes of the or¬ 
der of the symbol size. We therefore can safely conclude 
that our finite size results support a finite SHC at low 
densities even for strongly disordered systems. 

Naturally, the system gets more stationary with disor¬ 
der, as it is shown in panel (c) of Fig.|^ where a small 
residual value of 7 for the clean striped system (black 
curve) is suppressed for all /r’s. Again, for /r < —4t this 
does not imply the vanishing of cr®^, as would be the 
case for homogeneous Rashba SOC. 

— Discussion and conclusions — The above analy¬ 
sis clearly demonstrates that a system with modulated 
Rashba SOC can sustain a finite SHE in stationary con¬ 
ditions. This occurs for a limited density range, when the 
chemical potential falls in a region where the states are 
strongly affected by the modulated a{x) and are almost 
localized in the bottom of the modulating potential (in 
the direction of the modulation; the states are extended 
in the perpendicular direction). Therefore the response 
of the charge current to the electric field along the 
modulation direction is strongly suppressed which can 
lead to large spin Hall angles for the striped 

system. It is important to note that this effect is due to 
the modulation of the Rashba SOC and cannot arise in 
a “conventional” charge density wave. In fact, for con¬ 
stant a and independent of the electronic structure Eq. 
([^ predicts the vanishing of the SHE under stationary 
conditions. 

The implementation of this analysis in a real system is 
for sure a challenging task for several reasons. First of all 
the top gating structuring has to be sharp enough to pro¬ 
duce a sufficiently sharp spatial modulation of the 2DEG 
below the LAO layer (which is at least 20nm thick): if 
the modulation of the SOC is not sharp enough on the 
Lso scale, the 2DEG would feel a nearly uniform a and 

















5 


the SHE is expected to vanish. One should also con¬ 
sider that our analysis is based on a simple one-band 
model, while the SOC in the 2DEG in the LAO/STO 
involves several bands [iiiniisi]- Of course, the basic 
ideas of this work could be tested and hopefully imple¬ 
mented in other, perhaps simpler, systems involving het¬ 
erostructures of semiconductors with modulated Rashba 
SOC which have been already discussed in the literature 
in different contexts |55H55] . 
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edges support from the Deutsche Forschungsgemein- 
schaft. M.G. and S.C. acknowledge financial support 
from Sapienza University of Rome, project Awards No. 
C26H13KZS9. 


APPENDIX 


Rashba systems in the diffusive limit 

To gain insight on the numerical results discussed in 
the main text, we discuss here a continuum Rashba model 
in the diffusive limit. The corresponding diffusion equa¬ 
tions can be derived from a microscopic model by using 
e.g. the Keldysh technique. Such a derivation is, for in¬ 
stance given in Ref. IMl For the following discussion, 
however, one does not need to know such a microscopic 
derivation in detail. One main advantage of the diffusion 
equation description is that it contains all the important 
aspects of the Rashba model and allows an almost an¬ 
alytic treatment, which helps in elucidating the physics. 
We first provide the diffusion equation description for the 
uniform case. This is very standard and a recent discus¬ 
sion can be found in Ref. [3^ Then we describe the 
single-interface problem as the simplest realization of a 
non-uniform Rashba system with two regions with dif¬ 
ferent Rashba SOG. A subsequent subsection reports the 
two-interface problem. 


The uniform case 

In the presence of the Rashba spin-orbit coupling, the 
spin polarization along the y direction obeys the follow¬ 
ing equation in the diffusive regime (see Ref. [36] for a 
derivation) 


-dtsy+ Ddlsy =—{sy-so), (s) 

Top 

where D = v\tI2 is the standard diffusion coefficient. 
T^p = {2ma)^D is the inverse Dyakonov-Perel spin re¬ 


laxation time. So = —^NoarEx is the non-equilibrium 
spin polarization induced by the electric field applied 
along the x direction. Such a non-equilibrium polariza¬ 
tion is sometimes called the Edelstein effect or the spin- 
galvanic effect. Here Aq = m/2Trti^ is the 2D density 
of states and, in the following, we take units such that 
h = 1. We consider only the dependence on the x di¬ 
rection in the diffusion equation to make contact with 
the numerical calculation. In stationary and uniform cir¬ 
cumstances, we must have s*' = Sq. This leads to the 
vanishing of the spin current Jy, as it is well known in 
the Rashba model (see for instance Ref. [35|l . To see this, 
consider that the spin current is given by two terms: a 
” drift-like” Hall term and a ” diffusion-like” one. The 
drift term corresponds to the calculation of Ref. [T31 i.e. 
it is just the Drude formula for the spin Hall conductiv¬ 
ity. As for the ordinary Hall conductivity, it is non-zero 
and finite even in the absence of disorder. The diffusion 
term is usually expected to vanish in uniform situations. 
However as derived in Ref. |33| and used in Ref. [331 the 
Rashba interaction can be described in terms of a SU(2) 
vector potential, which then introduces covariant deriva¬ 
tives. The latter are defined by 





(9) 


where Af is the SU(2) vector potential. In the Rashba 
case. Ay = —A^ = 2ma. The key observation is that the 
diffusion term related to the covariant derivative of the 
spin density is present even in the uniform case. Specifi¬ 
cally the spin Hall current we are interested in reads 


TZ _ p I ^ ^ 

Jy-ao E- + 


( 10 ) 


where Lso = (2ma)“^ is the spin-orbit length. Notice 
the relation L'^q = Dts, which amounts to say that the 
Dyakonov-Perel relaxation time is the time over which 
electrons diffuse over a spin-orbit length, ctq^ corre¬ 
sponds to the expression given in Ref. [T31 i.e. the Drude 
formula evaluated without vertex corrections, of the spin 
Hall conductivity 


a 


sH 

0 


e 2r 
Stt Top ' 


( 11 ) 


In the diffusive regime appT 1, one has top ^ t, 
whereas in the ballistic one appT ^ 1 top t- Notice 
that, by using Eq. 0 and the expression for sq) the spin 
current (10) can also be written as 


Jy = 7^-isy-So). 
^ 2ma To 


( 12 ) 


By identifying the y-polarized spin current flowing along 
the X direction 


jy = -Ddxsy 


(13) 
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the diffusion equation becomes 


After solving the system for Ss- and (5s+,we get, 


dtsy+ d,jy = -2maJ^ = --^J^ (14) 

Lso 

which is the continuity-like equation for showing how 
the torque term associated to is expressed in terms of 


Single-interface problem 

In this Section we consider a static time-independent 
situation. The idea is to analyze a single-interface prob¬ 
lem in order to understand the supercell numerical calcu¬ 
lation. We then assume the following expression for the 
Rashba coefficient 

a{x) = 0{x)a++ 9{—x)a-. (15) 

Clearly at ±oo, one recovers the uniform case with com¬ 
plete cancellation of the spin current for the Rashba 
model with couplings a±. The strategy is to solve the 
diffusion equation in the two regions and connect them 
via the appropriate boundary conditions at a; = 0. We 
seek a solution of the form 

sy{x) = 9{x) (^so,+ -I- (5s+e““/'^+^ 

+ 9{-x) (^Sq-+ . (16) 

In the above so,± are the asymptotic values of the y- 
spin polarization at ±oo. Also L± are the corresponding 
spin-orbit lengths in the two regions. The z-polarized 
spin current in the two regions reads 

^ (so.± + 

= (17) 


6s- = Asi 


D_ 1 

' A- D+L-^ + D_LZ^ 


5s+ = —Asi 


D. 


L+ D+L+^ 


D_L 


-1 • 


( 20 ) 

( 21 ) 


From this it is clear that the spin Hall current averaged 
over the spin-orbit length 


1 

ZI 


dx J,! 


1 


dx Jy = 0, 


( 22 ) 


which is the analog of Eq. (14). Furthermore the left-hand 


side correspond to the quantity 7 in the main text. On 
the other hand we have that the total z-polarized spin 
current is different from zero 


dx Jy _ + 


dx Jy + = 


(23) 


2eNoD+D_mia+- a_)^T _ .h r p 

’’-^x — ^ ^ef f 


D.L 




D_L- 


where Te// is an effective length determined in terms of 
L± and D±. The spin Hall current is localized at the 
interface within a distance of order L±. 

The above calculation of a single interface suggests 
that the calculation for the periodic modulation de¬ 
scribed in the main text can be analyzed in terms of a 
series of interfaces. The spin current flows in the interface 
regions. However, by making the interface separation of 
the order of the spin relaxation length, one may have a 
spin current finite everywhere. 


Two-interface problem 

Here we consider the problem with two interfaces with 
the model given by 


where in the last step we used the fact that the constant 
terms cancel in each region. D± may differ in the two 
region because via the Fermi velocity they depend on 
the electron density. Eq.([T7|) shows that in the interface 
region there can be a spin current different form zero. To 
evaluate it we need to know the two values i5s±. To this 
end we use the continuity of the spin density and of 
the spin current J| at the interface. Continuity of the 
spin density gives 


5s- — 5s+ = Asq, (18) 

where Asq = so,+ ~ so,- = —eNQT{a+ — a-)Ex. Conti¬ 
nuity of the spin current instead gives 

+ ^6s+ = 0. (19) 


a(x) = a+9{l — |x|) -I- a_d(|x| — 1). (24) 

The solution for the spin density (x) is of the form 


sy{x)= 9 {x — l) so-+Ssiie 

+ 9{-x-l) L-+<5sLe("+')/^- 


+ 9{l-\x\) 


So,-I- ~ Aso 


Ss 


sinh((Z — x)/L_|.) 
sinh( 2 l/L-|.) 


cosh(x/L+) 

cosh(Z/L+) 

sinh((^-Hx)/L+) 
^ sinh( 2 Z/L+) 


where the continuity of s*' has already been implemented. 
The two constants Ssj^ and Ss^ must be determined by 
imposing the conservation of the longitudinal y-polarized 
spin current jy as in the single-interface problem. 
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After going through steps as for the single-interface we 


get 


/ OO 

da; Jy{x) = (25) 

-OO 

4eA^0-O+-D_m(Q;_|_ — Q;_)^r tanh(Z/L_(_) 


tanh(^/L+) -|- D-L_ 


-1 




The above equation has a simple interpretation in the 
limit I ^ OO, when the two interfaces are far apart. In 
this limit the two interfaces are independent. The total 
spin current is the sum of the spin currents flowing at 
the two interfaces. Eq.(25) reduces indeed to twice the 
contribution (23) of a single interface. On the other hand 
when I ~ L-|. the two interfaces interact and the spin Hall 
current in non zero everywhere. 


Continuity equations for the Rashba lattice model 

We provide here a discussion to clarify the role of an in¬ 
homogeneous Rashba SOC to allow for a SHE in station¬ 
ary conditions. The electron spin S obeys the equation 
of motion 

f = -.|S,Hl (26) 

and for the following it is convenient to separate the com¬ 
mutator into a term related to the divergence of spin 
currents and a ’rest’ G 


[S, H] = -idivJ + iG. 


(27) 


As a result Eq. (26) can be interpreted in terms of a 
continuity equation 

dS°‘ 
dt 


G“ = divr 


(28) 


where G acts as ’source’ term which in general is finite 
due to the non-conservation of spin. 

From evaluation of the commutators one finds the re¬ 
lation 


i.e. a torque for the y-component of the spin is associated 
with a z-spin current when the Rashba SOC a 0. Upon 
combining Eq. (29) with the y-component of Eq. (28) 
one finds 

dS^ 

|- [divJ^]i Qfiy-i-y Jj j_|_y -|- ai—y^iJi_y j^ = 0 (30) 
and in particular for a system with periodic boundaries 

f dS^ I 

^ ' I ^ ^i,i+yJi,i+y A ^i-y,i'^i-y,i j ~ ^ 


since the total divergence of any current has to vanish. 


For a homogeneous SOC coupling Eq. (31) implies that 


the total z-spin current has to vanish under stationary 
conditions. 

We exemplify the situation for a inhomogeneous 
Rashba SOC which varies along the x-direction as 


1 

Cti,i+x — 2 


ao + ai + (oo - ai ) sgn(sin ) 


(32) 

(33) 


i.e. one has stripes of width L with Rashba SOC oq alter¬ 
nating with L-wide stripes having coupling oi (cf. Fig. 
[^. The inhomogeneity in ai^i+x leads to a concomitant 
charge modulation which is shown in Fig. [^for the case 
L = 4. Clearly the charge accumulates in regions with 
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FIG. 5. Stripe-like modulation of the Rashba SOC along the 
x-direction. Stripes of width L and coupling oq alternate with 
stripes of the same width and coupling ai as indicated by thin 
and thick bonds, respectively. 


large Rashba SOC leading to a CDW profile. We can 
therefore view this model as an ’effective’ model for a 
density dependent Rashba SOC. 
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FIG. 6. Modulation of the coupling constant ai,i+rc (red) and 
charge density along the x-direction for the L = 4 stripe-like 
Rashba SOC. Particle concentration: n = 0.07. 


Fig. 0 shows the currents flowing in the (stationary) 
ground state. Thus from Eq. (28) the torques Gf are 
completely determined by the divergence of the spin 
currents which are shown by squares in the top panel 
of Fig. As a consequence of Eq. (29) a large z-spin 












































current is flowing on sites where also the y-torque is large 
in contrast to a homogeneous system where = 0. In 
the ground state the total y-torque vanishes so 

that from Eq. ( [2^ one obtains 
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FIG. 7. Spin currents (arrows) and torques (squares) in the 
ground state for a system with alternating L = 4 stripes with 
Rashba couplings ao = 0.2 and ai = 0.8, respectively. Top 
panel: y-components; Bottom panel: z-components. Particle 
concentration: n — 0.07. 


The total number of aoq-stripes for a x Ny lattice 
is Ustr — Nx/{2L). Denote with J§ ^ the total z-spin 
current flowing along the bonds of the ao,i-stripes. Then 
we can rewrite Eq. (34) as 


0 = aoJo"+aiJf ^ = (35) 

ai 

and the total z-spin current of the system is thus given 


by 


= Ustr {JS + JD = nstrJS ( 1 - - 


(36) 


Now we can draw the following conclusions: The mod¬ 
ulated Rashba SOC causes local torques G\ which are 
related to local flows of z-spin currents. Provided that 
the total y-torque vanishes the system thus exhibits a net 
flow of Jfgj for oq 7 b oi in the ground state. 

The same analysis can now be applied in the presence 
of an electric field = —dtA^ which couples to the sys¬ 
tem via the charge current, i.e. H' = —e t) 

and Obviously Eq. (30) holds also in the result¬ 


ing non-equilibrium situation which we eva luate in linear 
response. Each operator Oi in Eq. (30) reacts to the 


electric field according to Oi{uj) = ieKij{uj)Ej{ijj )/uj and 
the correlation function is given by 


A,,(a;) = -- / df0(t-t')e“(‘-‘')( 

J —00 






FIG. 8. Top panel: Modulation of the coupling constant 
Oii^i+x (red) and induced J^-current (black, circles) for a cut 
along the x-direction and L = 4 stripe-like SOG. The blue 
line (squares) shows the temporal change of the y-spin com¬ 
ponent dSy{ix)/dt. Bottom panel: Induced spin current 
(black, open circles) and corresponding divergence (black, 
full circles). The blue line (squares) shows the response of 
yij; = 2ai^i+yJi i^y along the same cut. Particle concentra¬ 
tion: n = 0.07. 

The top panel of Fig. shows the induced spin 
current together with the induced temporal change of 
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dSy{ix)/dt along a cut in x-direction of a L = 4 striped 
system. The dominant contribution to comes from the 
boundary regions between small and large a stripes which 
gives rise to a finite spin Hall conductivity. Moreover, in 
contrast to the homogeneous case, where the induced 
current and dSy{ix) / dt are equal 


dS^ 2ma 


(37) 


(although with opposite sign), we now observe a much 
more stationary behavior. The reason can be deduced 
from the bottom panel which reports the x-dependence 
of the induced spin current along with its diver¬ 
gence, i.e. [divjy]i = Jfx.ix+ii^x) - J!x-u^x■ It turns 
out that the spatial behavior of the contribution ^ix = 
‘^Oii,i+yJii+y{ix) is similar to [divjy]i but opposite in sign 
which from Eq. (30) is responsible for the small value of 
dSy {%x ) /dt. 


From the above considerations we see that in a ho¬ 
mogeneous Rashba system a finite spin Hall conductivity 
necessarily implies a non-stationary situation with the lo¬ 
cal accumulation of spin density. On the other hand 
an inhomogeneous Rashba coupling partially shifts this 
time dependence to a finite divergence of the spin 
currents which is stationary due to non-conservation of 
spin. 
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